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B. Z. Kopeliovich, Ivan Schmidt and M. Siddikov 

Departamento de Fisica, y Centra Cientifico - Tecnologico de Valparaiso, 
Universidad Tecnica Federico Santa Maria, Casilla 110-V, Valparaiso, Chile 

We study the electromagnetic Bethe-Heitler type contribution to neutrino-induced deeply virtual 
meson production (j/DVMP). Such C'(Q:em)-corrections decrease with in the Bjorken regime less 
steeply than the standard i^DVMP handbag contribution. Therefore, they are relatively enhanced at 
high . The Bethe-Heitler terms give rise to an angular correlation between the lepton and hadron 
. scattering planes with harmonics sensitive to the real and imaginary parts of the DVMP amplitude. 

These corrections constitute a few percent effect in the kinematics of the forthcoming Minerva ex- 
periment at Fermilab and should be taken into account in precision tests of GPD parametrizations. 
For virtualities ~ 100 GeV^ these corrections become on a par with DVMP handbag contribu- 
tions. A computational code, which can be used for the evaluation of these corrections employing 
various GPD models is provided. 
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I. INTRODUCTION 



CD ■ 

[ Generalized parton distributions (GPD) allow evaluation of cross-sections for a wide class of processes, where the 
. coUinear factorization is applicable P, [3| ■ The main source of experimental information on GPDs has been so far the 
' electron(positron)-proton measurements performed at JLAB and HERA, in particular deeply virtual Compton scattering 
J> (DVCS) and deeply virtual meson production (DVMP) The 12 GeV upgrade at Jefferson lab will open new 

■"^ opportunities for further improvement of our knowledge of the GPDs p^ . 

However, the practical realization of this program suffers from large uncertainties. For instance, the results at moderately 
' high can be affected by poorly known higher-twist components of GPDs and distribution amplitudes (DA) of the 
produced mesons fl7l - l20j . 

Neutrino experiments provide a powerful tool for consistency checks for the extraction of GPD from JLAB data, 
' especially of their flavor structure. The study of various processes in the Bjorken regime may be done with the high- 
intensity NuMI beam at Fermilab, which will switch soon to the so-called middle-energy (ME) regime with an average 
neutrino energy of about 6 GeV, and potentially may reach energies up to 20 GeV, without essential loss of luminosity. 
In this setup the Minerva experiment [21] should be able to probe the quark flavor structure of the targets. Even higher 
luminosities in multi-GcV regime can be achieved at the planned Muon Collider/Neutrino Factory |22l424| . 

Certain information on the GPD flavor structure can be extracted from comparison of analogous processes in neutrino- 
. and electro-induced processes employing the difference of flavor structures of electromagnetic and weak neutral currents. 
An example is the weak DVCS [2^, which alone, however, is not sufficient to constrain the flavor structure. 

Recently we discussed the possibility of GPD extraction from deeply virtual neutrino-production of the pseudo-Goldstone 
mesons (tt, K, if) [26j . The i^DVMP measurements with neutrino and antineutrino beams are complementary to the 
electromagnetic DVMP. The octet of pseudo-Goldstone bosons, originating from the chiral symmetry breaking, acts in 
the axial current as a natural probe for the flavor content. Due to the V — A structure of the charged current, in z/DVMP 
one can access simultaneously the unpolarized GPDs, H and E, and the helicity flip GPDs, H and E. We expect the 
contributions of the GPDs Ht, Et, Ht and i?T, which are controlled by the poorly known twist-3 pion DA 0p, to be 
negligible. Besides, important information on the flavor structure can be obtained by studying the transitional GPDs in 
the processes with nucleon to hyperon transitions. As was discussed in [27], assuming SU{3) flavor symmetry one can 
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relate these GPDs to the ordinary diagonal ones in the proton. 

In this paper we study the Bethe-Heitler (BH) type radiative corrections to the diffractive neutrino-production of charged 
pseudo-Golstone mesons with the target remained intact, related to meson emission from lepton line with subsequent 
electromagnetic interaction with the target. Although such processes are formally suppressed by aem-, at large they 
fall off less steeply than the DVMP cross-section. While at virtualities < 10 GeV^, relevant for modern neutrino 
experiments, this is a few percent correction, already at ^ 100 GcV^ it becomes comparable with the DVMP cross 
section. Such corrections are relevant only in case of the i^DVMP: for the electron-induced DVMP ep -> epM they are 
suppressed by factor ^ (^GpQ^) , where Gp is the Fermi coupling, and are negligible unless we consider extremely high 
« Myy. Existence of such diagrams opens a possibility to probe separately the real and imaginary parts of the DVMP 
amplitude (not only the total cross-sections), in close analogy to DVCS studies [Tsl. [28l. [29j . 

The paper is organized as follows. In Section |TT] we present the results for the DVMP and BH contributions to the 
experimentally measurable total i/DVMP cross-section (technical details of evaluation may be found in Appendix[X|. Also, 
at the end of Section |TT] we construct two asymmetries which have particular sensitivity to the real and imaginary parts 
of the DVMP amplitude. In Section Hill for the sake of completeness we discuss the features of the GPD parametrization 
used in calculations. In Section ITVl we present the numerical results and make conclusions. 

II. CROSS-SECTION OF THE !/DVMP AND BH PROCESSES 

Exclusive neutrino-production of pious is presented by the diagram (a) in the Figure [T] Production of pious by the 
vector current was studied in ,31,] , and was recently extended to neutrino interactions in |2gJ . 




FIG. 1: Diagrams contributing to the neutrinoproduction of mesons, (a) DVMP process (b,c) BH contributions. 

As was shown in [sO, IsH , at large where the collinear factorization is applicable, the amplitude of this process is 
suppressed due to hard gluon exchange in the coefficient function, and is small. This raises a question, how important 
could be the O {aem) corrections ? While a systematic study of radiative corrections is beyond the scope of the present 
paper, we would like to consider in detail the contributions which decrease with less rapidly than the diagram (a). In 
the leading order in and aw there are two such diagrams (b,c) shown in the Figure[T] These diagrams are enhanced by 
a factor ^ /t, where —Q^ = and t are the lepton and proton 4-momcntum transfers squared, respectively. This factor 
is large in the Bjorken regime of i ^ Q^, and despite the formal suppression by aem, the diagrams (b,c) are numerically 
comparable to the diagram (a). Notice that such corrections are sizable only in the case of neutrino-production of mesons. 
In case of electroproduction the corresponding corrections are suppressed by ~ G\Q'^, where Gf is the Fermi constant, 
and are negligible, unless we go to extremely large ~ M^, . In what follows we evaluate the contribution of each 
diagram in Figure [TJ 

The cross-section of pion production can be presented as a sum of contributions of DVMP (diagram (a)), BH mechanism 
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(diagrams (b,c)) and their interference, 



dt dlnxBj dQ'^dfj) dt d Inxs j dQ^dcj) dt dlnxgj dQ^dcj) dt dlnxgj dQ^dcj)' 

where xbj = Q^/(2P • q), and (f> is the angle between the lepton and hadron scattering planes. 
Evaluation of the diagram (a) is straightforward and yields [2^ 



(1) 



^i^(DVMP) ^ Glflixlj (l 



2 2 
1^ 



dtdlnXB, dQU<^ g^^,^, ^ ^,^^^^^2 ^ 4^1^^ 3/2 



^ TmI', (2) 



where y is the fractional loss of lepton energy defined ds y — P ■ q/P ■ k — vjEy. Notice that the DVMP cross section 
turns out to be independent of the angle between lepton and hadron planes. This happens because the momentum q 
does not have transverse components in the Bjorken reference frame. 

For unpolarized target, the matrix element squared \Pm\ in Eqn. ([2]) can be simplified as. 



2 ^ 647r^ ^ 



4 (1 - XB, ) [Uu-Wu + UmUI^ - ^^mS;, (3) 



x%j {Um£*m + £mT-L\i + Um£*m + SmUXi^ - [x%j + (2 - xujf j £m£ 
where we introduced a shorthand notation, 



'-N 



M 



dz^^'- f\z^, (4) 

and the script letters H, £, H, £ signify convolution of the GPDs H, E, H, E with corresponding coefficient functions 
given in Table HI 



TABLE I: List of the DVMP amplitudes Hm, £m, Hm, £m for different final states. For a neutron target, in the r.h.s. we flipped 
Hu/n ~^ H^i/p, -ffd/n — > H^/p, so all the GPDs are given for a proton target. To get £, T-L, £ one should replace H with E, H and E 
respectively. Vij are the CKM matrix elements. c± is a shorthand notation c±(x,^) = l/(x±^=p iO) for the leading order coeflicient 
function. The NLO corrections to the coefficient functions may be found in [32 . |33 | . For the sake of brevity we do not show the 
arguments (a;, ^, t, Q) for all GPDs and omitted the integral over the quark light-cone fraction J dx everywhere. 



Process 


type 


Km 


Process 


type 


■Hm 


Vp fl~TT'^p 


CG 


Vud (Hdc- +Huc+) 


u n H~'K^n 


GG 


Vud {H^c- + Hdc+) 


up — >■ n'^Tv^p 


CG 


Vud {H^c- +HdC+) 


vn ^ H'^TT^n 


GG 


Vud (Hdc- +Huc+) 



up ^ ^ K^p 


GG 


Vus (C+Hu+C-Hs) 


un ^ H K^n 


GG 


Vus {c+Hd + c-Hs) 


up fi^K~p 


GG 


Vus {HuC- +HsC+) 


un ^ jj.'^K~n 


GG 


Vus {HdC- +HsC+) 



Comparing different elements in Table U one gets relations, 

j^DVMP _ 
""pp->/j+;T-p ^ ' 



A^DVMP _ j^DVMP [tr\ 
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which are just a manifestation of the isospin symmetry. As shown below, these relations are broken by the BH corrections. 

In the leading order in both BH diagrams, Fig. [T] (b) and (c), acquire dominant contribution from longitudinally 
polarized photons. However, as discussed below, certain angular harmonics are suppressed by ~ Aj^/Q and get similar 
contributions from transverse and longitudinal photons. For this reason, for the form factors we include both the longi- 
tudinal and transverse components and evaluate the BH diagrams exactly, and only after that we make an expansion in 
. The diagram (b) contains the matrix element 



1 



(7) 



where V^^(x) and A°(a:) are the vector and axial vector isovector currents. The correlator ([7|) can be evaluated in pQCD in 
the collinear approximation, because is large, and we assume that the dominant contribution comes from the leading 
twist-2 pion DA. Notice that the amplitude A^l.v cannot be interpreted as the pion form factor because (i) the virtuality 
of W is large, (ii) the insertion of the pion state between A° and leads to ^ q^, which gives zero acting on the 
transverse on-shell lepton current. The details of the calculations are moved to Appendix El and we present here the final 
result, which reads. 



where 



/ . o o ^ 5/2 



(8) 



riBH _ p 



BH 
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- ( (^4(0_i-l 
5t 



t 



1) 



(9) 



2m\ 



,2 



-1 + 5)) v^~A - 1)2 + ^ (40^1 - li^^i + 10)) y 



2m% 



1 - 2)' + - 1)') xl + {2y - 1) ^ (0-1 - 1) (-(/--i + y (20_i - 1) + 2) xg 
- 4 (1 - 2yf - 1)^) F^it) + 2F,{t)F,it)xl (xj, (y + if - ^ (y + if 
V ((8</'-i - 240-1 + 17) - 2 (802^1 - 240_i + 19) y + lO^^^i - 360_i + 35) ) 



7^ 



I 

+2 (^6(f,\ - 2O0_i + ^ (20_i -1) + I7^y- 94>\ + ^ (5 - 40_i) + 340_ 
- {2y - I) 1^-^^ (0_i - 1) (-</)_! + y (20_i - 1) + 2) + (1 - 2yf [j^^ 



34' 



6-1-1)' 
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BH 



Km 



N 



9Q2 

2(0-1-1 
it 



4(3(-4y + 3(2/-2)0_i+9)a; 



(10) 



2^ 



2y 



t(f>-i 



m 



N 



3t 



+ 30-1-6) +18)x| 



6t 



— (0_i - 1) (-60_i + y (40_i -3) + 3)xB- {2y - 3) — (0_i - 1)^ Fl{t) 



'TV 



+ AF^{t)F2{t)x'B ( 3 (y - 3) a:^ - 12 (y - 3) - 8^ (4y (0_i - 3) - 50_i + 18) (0_i - 1) 



FKt) (-6(y-3)x| 



+ (-6 {y - 3) 4 + 12 (-2y + 3 (y - 2) 0_i + 3) + 8 (2y (0_i - 6) - 0_i + 18) (0_i - 1)) : 



+24 



/ t 



V4m2, 

' ' Q4 



(XB (XB - 20_i + 2) + 2 I 



1)){XB (y-3)-2(2y-3) (0-i - 1)) 



4 if 2 



1 



[(5X5 



4m^ 
t 



^1+4) (0_i 
4m|r 



l)F2(t) + 2x|Fi(t)^^2(t) 



1-1) 



'AT 



Here, following Refs. [H, [H, [2^ we introduced the notations, 



(11) 



1-y- 



1 - 



I - XB 



o\ — 



+ + 



4xi3(l-a;s) + e^i-<,: 



4(1-xb) 



(12) 



(13) 



Notice that in the leading order in the contributions of diagrams (b) and (c) to the terms Cq^ and Cf^ exactly cancel 
each other, so that the coefficients C^^ and C^^ acquire an extra suppression factor m\/Q'^ . This can be understood as 
a screening of the opposite charges of the pion and muon, when both of them move in the forward direction in the limit of 
massless leptons. As we can see, the BH cross-section is symmetric under the 4> — >■ —0 transformation. For asymptotically 
large Q^^ the C^^ harmonic is suppressed by A_l/Q, whereas C^^ ~ C^^, so the distribution is symmetric relative to 
the replacement tt — 0. Note also that K ^ and vanishes when A_l 0. 

The interference term in ([TJ has a form (see details in Appendix |^ 
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fiG% XBa,mas(l)-i {Ct + Cf* cos + S^' sin , 



dtd\nx Bj d.Q'^ 



36ttHQ^1- ^) 1 



5/2 



(14) 



where 



nint 
^0 



'AT 



V 4m^ 



X ( {2^_^-i)x%- — {^^,-l){l + XB) ]+0[-^,— 



if 



Ci"* = y [^1 (2 7^e£: {y - 3) x| + 7^eH (4 (xb - 2) {2y - 3) - 4 ((xs - 4) y + 6))) 



+ F2 2 7^e-H(y- 3)a;^ +7^e£ 2x| (y - 3) - (0:5 - 2 



4m^ 



(4(2y-3) (0_i-l)-2xB (y-3)) 



t 

K{2-xb) 
6 



[i^i {-2Im£{y + 1) x| - 2Tm?^ (xs - 2y - 2(/._i + 4) - 4 (2y - 1) - 1))) 



F2 {-2ImH{y + l)xl 



2Im£ (y + 1) 1 



' ' Q2 



4m2r 



2m 



~ (-2(/._i2/ + + -2)xB + {2y - 1) ^ 



-1) 



AT 



'AT 



(15) 



(16) 



(17) 



As one can see from (jl4p , the angular dependence of the interference term has a sin (j) term which is absent both in the 
BH and DVMP taken alone. This term stems from the interference of the vector and axial vector current in lepton part 
of the diagram. It has the same sign for neutrino and antineutrino beams (in contrast, all the other terms in (1141) change 
sign). Such terms are absent in case of DVMP and BH since the interference is asymmetric w.r.t. polarization vectors of 
the emitted boson. Due to presence of sin cf) harmonics, the antisymmetrized cross-section directly probes the imaginary 
part of the DVMP amplitude as 



dtd\x\XBj dQ'^d(f) dt dlnxBj dQ^dcj) dt dlnxBj dQ'^dcj) 

Another way to access the interference term is based on the isospin symmetry for the pure DVMP cross-sections (j5l6p . 
Since the BH correction nearly vanishes on a neutron target, and that BH cross-sections may be easily calculated, one 
may directly probe the interference term and extract the real and imaginary parts of the DVMP amplitude as, 



d'^T^"*^ + dV- + d^a + dVJ'"'> + 



dtdlnxBj dQ'^dcj) dtdlnxBj dQ'^dcj) dtdlnxBj dQ'^d<p dtdlnxBj dQ'^dt 



dtdliiXBj dQ'^dcj) dtdlnxBj dQ'^d(j) dtdlnxBj dQ'^dcj) dtdlnxBj dQ'^ 



(19) 



(20) 
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III. GPD AND DA PARAMETRIZATIONS 



As was mentioned in the introduction, an essential uncertainty in the calculations of DVMP originate from the poorly 
known DAs of the produced mesons. Only the DAs of pions and ry-mesons have been tested experimentally, and even in 
this case the situation remains rather controversial. The early experiments CELLO and CLEO [sj], which studied the 
small-Q^ behavior of the form factor FM-yy, found it to be consistent with the asymptotic form, <f>as{z) = 6 z{l — z). Later 

the BABAR collaboration [s^ found a steep rise with of the form factor IFt^^^ (Q^) |^ ^^'^ large-Q^ regime. This 
observation gave rise to speculations that the pion DA might have a z-dependence quite different from the asymptotic 
form ^3^1 (see also the recent review by Brodsky et. al. in [s^, H^). However, the most recent mesurements by the BELLE 
collaboration did not confirm the rapid growth with observed in the BABAR experiment. As was found in ^40. l4l| 
based on the fits to BELLE, CLEO and CELLO data, the Gegenbauer expansion coefficients of the pion DA (I>2;tt{z) 
are small and give at most 10% correction for the minus-first moment defined in Although there are no direct 
measurements of the kaon DAs, it is expected that their deviations from the pion DAs are parametrically suppressed by 
the quark mass ms/GeV . Numerically this corresponds to a 10-20% deviation. 

For this reason in what follows we assume all the Goldstone DAs to have the asymptotic form, 

'/'2;{7r,A:,,,}(^) ~ (t>as{z) = 6z(l - z) . (21) 

For the decay constants we use the standard values « 93 MeV, fK~ 113 MeV, and fr, ~ Ik- 

For GPDs more than a dozen different parametrizations have been proposed so far fa. l4_ 42 - isj . While we neither 
endorse nor refute any of them, for the sake of concreteness we select the parametrization [42l |49|. 50| . which succeeded to 
describe HERA [ST] and JLAB [4^, |4^, [H^l data on electro- and photoproduction of different mesons, so it might provide 
a reasonable description of vDVMP. This parametrization is based on the Radyushkin's double distribution ansatz. It 
assumes additivity of the valence and sea parts of the GPDs, 



H{x,^,t) ^ H,aiix,^,t) + H,,aix,^,t), (22) 



which are defined as 



J\a\ + \I3\<1 ^['^-\P\) 



Hl^ = I dl3daS ~ X + " ' ^ g.aK/?)^^'''""' (23) 

|a| + |/J|<l 

HL. = f dpdaSiP-x + aO '-^^^^^^^^^^ (24) 

J\a\ + \I3\<1 8(1-1^1)5 

and Qyai and g^ea sue the ordinary valence and sea components of PDFs. The coefficients 6^, a^, as well as the parametriza- 
tion of the input PDFs q(x), Aq{x) and pseudo-PDFs e{x), e{x) (corresponding to the forward limit of the GPDs E, E), 
are discussed in [i^, 11^ [Eo]. The unpolarized PDFs q{x) within the range of < 40 GeV^ roughly coincide with the 
CTEQ PDFs. Notice that in this model the sea is flavor symmetric for asymptotically large Q^, 



where 



ff.". = «i. = ''(0')«i.. (25) 

The equality of the sea components for u and d quarks in (I25p should be considered as a rough approximation, since in 
the forward limit d ^ u was firmly established by the E866/NuSea experiment [52|. For this reason the predictions made 
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with this parametrization of GPDs for the p ^ n transitions in the region xsj G (0.1. ..0.3) might shghtly underestimate 
the data. 

The Dirac and Pauh form factors Fi{t), F2{t) are extracted from GPDs in the standard way, 

Fr{t)=J2ej\xHg{x,U), (27) 

q •'-I 

Fr{t)=J2eJ\xEgix,U)- (28) 

q •'-1 

IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we perform numerical analysis of the electromagnetic corrections to the processes hsted in Table U relying 
on the GPDs described in the previous section. At small the angular harmonics are small and the cross-section is 
dominated by the angular-independent DVMP contribution. Therefore, it is convenient to normalize all the coefficients 
to DVMP cross-section, 



dt d\nxB dQ^dcf) dtdlnxsdQ 




sismcj) . (29) 



Notice that in the limit aem — >■ 0, no BH corrections are possible, the coefficient cq = 1, and all other coefficients vanish. 

The results for the Q^-dependence of the relative BH corrections to the neutrino-DVMP cross section for pions and 
kaons are presented in Figure [2l 

We see that the isotropic part of the BH correction 1 — cg steeply rises in all channels from few percent or less at 
< lOGeV^ up to few tens of percent at ~ 100 GeV^. It behaves like 1 — Cg oc modulo logarithmic corrections. 
As a result, the cross section is reduced about twice compared to the DVMP contribution. The asymmetry si also rises 
with and reaches about 15 % at = 100 GeV^. Notice that some of the coefficients (e.g. Cq — 1, Ci) have nodes for 
TT'^ /K^ production, while they are absent for tt~ /K~. The reason is purely algebraic: as one can see from the Table U 
for 7r+/_ft'+ the large s-channel coefficient function c_ is convoluted with the small d/s -quark GPD, whereas the small 
M-channel coefficient function c+ is convoluted with the large u-quark GPD. This produces a node in the real part of the 
DVMP amplitude, because the real parts of the two contributions have opposite signs. Such a node is absent for negatively 
charged mesons, because the "large" w-quark GPD is convoluted with the "large" c_. The full DVMP cross section has no 
nodes, because it gets a large contribution from the imaginary part, which homogeneously depends on (the coefficient 
si, which probes the imaginary part, has no nodes). The difference between the Cabibbo suppressed and allowed processes 
comes from the sensitivity to different flavor combinations of GPDs in the corresponding DVMP amplitude. 

The terms Cp — 1 and Ci in Eqn. (I^^l are dominated by the interference of the DVMP and BH amplitudes, therefore 
they have different signs for 7r+ and tt" (and and K^). The term C2 gets contribution only from BH process, so it 
always has the same sign. The si-term does not change its sign under the C-conjugation in the lepton part, because it 
originates from the P-odd interference between the vector and axial vector currents. 

The results for Aj^-dependence of the BH corrections are depicted in Figure |3l 

It shows that there is a qualitative difference between cq and other angular harmonics. The coefficient cq reaches its 
maximum at A_l — due to the 1/t behavior of the BH cross section. In contrast, the angular harmonics ci, C2, si vanish 
at small Aj_ due the if-factors in front of them. As a consequence, the harmonics reach their maxima at Aj^ ^ 0.1 GeV. 

The results for the angular harmonic coefficients vs the elasticity parameter y = P ■ q/P ■ k = vjE^, are presented in 
Figure IH 

The coefficient C2 does not depend on y at all due to exact cancellation of the pre-factors 1 — y — y^e^ /A in and 
the pre- factor in the DVMP cross-section Eqn. ©. The harmonics cq and ci have a mild dependence on y, except in the 




region y ^ 1, where they blow up, because the DVMP cross-section ([2]) is suppressed there by the factor ^ 1 — y — 
whereas the harmonics in (j8ll4l) are suppressed at most as K oa y^l — y — y'^e^ /A. The harmonic si in accordance with 
Eqn. pT)) in the kinematics xb <C 1, \t\ <C mj^ is proportional to ^ (2y — \)Fi{t)ImT-L and has a node near y 0.5. 



V. SUMMARY 



In this paper we studied the electromagnetic Bethe-Heitler corrections to neutrino-induced deeply virtual meson pro- 
duction. We found these corrections to fall with less steeply compared with the z^DYMP cross section, so they tend 
to become a dominant mechanism in the Bjorken limit of — >■ cjo. Besides, they are enhanced at small-t due to the 
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FIG. 3: (color online) Ax-dependence of the BH correction to the j/DVMP process. 



t-channel Coulomb pole cx 1/t. Remarkably, these corrections generate an angular correlation between the lepton and 
hadron scattering planes. Similar to the BH corrections in DVCS, some angular harmonics are sensitive to the real or 
imaginary parts of the DVMP amplitude (see (|15m7p ). Notice that the appearance of such angular dependence was also 
predicted in (53i] , however there it appears due to interference of the longitudinally and transversely polarized charged 
bosons, which is a twist-three effect. In contrast, our result is a twist-two effect, and as one can see from Figure [5J it is 
not suppressed for asymptotically large Q^. 

Numerically, the BH corrections are subject to the interplay between the suppression factor ^ Ofem and the relative 
enhancement, which is as large as ~ Q^/t for some harmonics. In the kinematics of the Minerva experiment, the BH 
contribution for the proton target represents a few percent correction and thus is important for precision tests of the 
GPD parametrizations. At ~ 100 GeV^, which can be accessed in future neutrino experiments, these corrections 
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FIG. 4: (color online) y-dependence of the BH correction to the z^DVMP process. 



are expected to become on par with the DVMP contribution. For a neutron target, these corrections are two orders of 
magnitude smaller than for a proton and can be neglected up to very high Q^. Combining this fact with isospin symmetry 
of the DVMP amplitude, we construct combinations of the cross-sections Eqs. (|19l20p . which are sensitive only to the 
interference term. 

The electromagnetic corrections discussed in this paper are important only in neutrino- induced DVMP: in the case of 
electron-induced processes ep ^ epM the BH corrections are suppressed by the factor (Gj^Q*), and are negligibly small. 
We provide a computational code, which can be used for evaluation of the cross-sections relying on different GPD models. 
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Appendix A: Evaluation of the i^DVMP and BH cross-sections 



In this section we present some technical details of the evaluation of diagrams (a-c) in the Figur e IT] T he calculation of 
the diagram (a) in Figure [T] is rather straightforward and yields for the amplitude of the process 26l l30l l3ll | 



Ta 



STTi Gf M(fcM)gI(g)(l-75)t^(fc|.) 



V2 



Q 



(Al) 



where p.{k) and v{k) are the spinors of the final muon and initial neutrino; e{k) is the polarization vector of the photon; 
N{p), N{p) are the spinors of the initial/final state baryons; 4)m{z) is the DA of the produced meson; /m is the decay 
constant of the meson M] subscript index for each momentum in Eqn. (jAll) and in what follows shows to which particle 
it corresponds; Hf^N {p^) TN (pi) is a symbolic notation for summation of all the leading twist GPDs contributions 
(defined below); and H\,f are the convolutions of the GPDs Hr of the target with the proper coefficient function. Currently, 



the amplitude of the DVMP is known up to the NLO accuracy [32, [32] ■ Extension of the analysis of [30|, |31| to neutrinos 
is straightforward. In contrast to electro-production, due to the V — A structure, the amplitudes acquire contributions 
from both the unpolarized and helicity flip GPDs. 

Four GPDs, H, E, H and E contribute to this process in the leading twist. They are defined as 



"2^ 



dz 



= {H,{x,^,t)N{p2)j+N{pi) 
Afe 



-Eq{x,i,t)N{p2)ia+kN{pi) 



iP2) Ipq' (-|) l+"/5^q (0 A{pi) 



2m N 

= [Hgix,^,t)N{p2)j+j5N{pi) 



(A2) 



(A3) 



2niN 



Eg{x,C,t)Nip2)N{pi) 



where P = pi + p2, A — p2 ~ pi and ^ — A+/2P+ xbj/ (2 — xbj) (see e.g. [iH for the details of kinematics). In what 
follows we assume that the target A is either a proton or a neutron. Since in neutrino experiments the target cannot be 
polarized due to its large size, it makes no sense to discuss the transversity GPDs Ht, Et, Ht, Et- We also ignore the 
contributions of gluons in this paper because in the current neutrino experiments the region of small xsj ^ 1 but very 
high Q^, is hardly accessible, so the amplitude (jAl|) simplifies to 



T = 



Sttz Gf M (^m) ^L{q)v {ky) Us 
{l + QVM^) Q 



dz 



4>m{z) 



^£mN{p2)i5N{p^) 
ZniN 



(A4) 



nMN{p2)-/+N{pi) 



Afc 
2mN 



£mN {p2)i(y+kN {pi) 



In Table U the corresponding amplitudes are listed for each final state M . The DVMP part of the corresponding neutrino 
cross-section for charged currents is given by 

In the leading order in both BH diagrams Figure [1] (b,c) are dominated by longitudinally polarized photons. Nev- 
ertheless, as was mentioned in Section |lll we evaluate the BH contribution exactly, because various angular harmonics, 
suppressed by ~ A±/Q, get contribution from the transverse components, which is of the same order as the longitudinal 
result. Only after that we make expansion in 

The dipole scattering amplitude, which contributes to the diagram Figure [T] (b) , has the form 



All (g, A) ^j-J d^xe-*^- {0\{V^{x) - A^ix)) Jr(0)| tt^ {q - A)) 



(A5) 
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where V^{x) and A^j^i^x) are the vector and axial- vector isovector currents. Notice that the amplitude should not be 
interpreted as a pion form factor, because: (i) the virtuality is large; (ii) the insertion of the pion state between and 
J^™ leads to ^ (7^, which vanishes when is multiplied by an on-shell lepton current. 

We evaluated (|A5p in pQCD in the collinear approximation. This is justified in Bjorken kinematics by the high virtuality 
of the charged boson, so we assume that the dominant contribution comes from the leading twist-2 pion DA. The result 
reads, 

Al^ij ^) I (5m"^/o + + 91^"^) /i + (^M^i' + A^n^) /2 - ief^^isyuisq^gi - ie^,yp^ni3/\^g2) ; (A6) 

f ^ ^ {-2zz + (1 - 2zz) t/Q') n-q-{{l- 2zz) - 2zzt/Q') n ■ A ^ n f M , 

f, /X 2zz- il~2zz)t/Q^ 2 ,^(m\ t\ , 

f , , (1 -2zz) -2zzt/Q^ 2 , , ,^(m\ t\ , 



3 g • n 7 ' ' XB (z - zVQ^) (z - zt/Q2) 3^;^ VQ^ ' 

Here and are the positive and the negative direction light-cone vectors respectively. The plus-components of q and 
A have the form, 

(xB - t/m%) (1 - \/rT^) - XB {2xB + 1 - VTTT^ t/Q^ 
n-A — - — , ; {A13) 



where e = 2mNXB /Q- In what follows we encounter the combination 4 — /o — xsifi — /2), for which we need to 
make expansion up to 0{Q~^). While separately the series expansion coefficients for each factor fi have non-integrable 
singularities ^ z^^z"^, which signal a sensitivity to the transverse degrees of freedom and presence of the non-analytic 
terms ~ \a(Q'^ /\t\)/Q'^, in the above-mentioned combination, these terms cancel each other resulting in 

/o + ih - /.) - 4 + ,^ Ji^ + -B) -^rn%xl _ 2t (1 + xb) - 6m%xl _ 

The amplitude of the axial current transition into an on-shell pion in the diagram Figure [T] (c) according to PC AC has 
the form, 

(0 I j;;'5(0)| 7r'^(g - A)) = iUV2q^ (A15) 
In order to simplify the calculation of the leptonic part of the diagram (c), we employ the chain of identities 

K (1 - 75) >^ {K) = (1 + 75) Ki' {K) = (1 + 75) [k^ - V (fc^) , (A16) 



14 

S [K - K) (1 + 75) (fc. -k)=S [K - K) [k - k) (1 - 75) (A17) 

= -1 + 75 + rnf_,S {k„ - A:^) w - (1 - 75) , 



Ai - e/i (k^) i {k^) S {ky - k^) k^ (1 - 75) v (k^) (A18) 
= -ep, (fc^) e (kj) V (ky) + O (m^) , 

where in (IA16P we make use of the fact that the initial state neutrino is on-shell, kyv(k) = 0. Actually, the simplification 
of (jA18|) is a manifestation of the Ward- Takahashi-Slavnov- Taylor identity for the charged lepton current, 

{flic. (1 - 75) i^) = -^^w^^ (1 - 75) + 7!^^^ " ^^'^ ' ""^^ ~ 



9\9v^9a -9v -9a 



2 cos e 



w 



/iZ (1-75)1^^0. (A19) 



This simplification explains why there is no harmonics in the denominator of the BH and interference terms, and it is 
valid in the limit of massless leptons. The diagrams (b, c) yield for the amplitude (sign corresponds to 7r+) 

Further evaluation of Bethe-Heitler ([5]) and interference terms requires some trivial but tedious Dirac algebra. 



which was done with FeynCalc 54 
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